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Abstract: Recently, suboptimality estimates for model predictive controllers (MPC) have been 
derived for the case without additional stabilizing endpoint constraints or a Lyapunov function 
type endpoint weight. The proposed methods yield a posteriori and a priori estimates of the 
degree of suboptimality with respect to the infinite horizon optimal control and can be evaluated 
at runtime of the MPC algorithm. Our aim is to design automatic adaptation strategies of the 
optimization horizon in order to guarantee stability and a predefined degree of suboptimality for 
the closed loop solution. Here, we present a stability proof for an arbitrary adaptation scheme 
and state a simple shortening and prolongation strategy which can be used for adapting the 
optimization horizon. 

Keywords: adaptive control, model predictive control, stability, suboptimality, sampled-data 
control, sampled-data systems 



1. INTRODUCTION 

Nowadays, model predictive controllers (MPC), sometimes 
also called receding horizon controllers (RHC), are used in 
a variety of industrial applications, cf. Badgwell and Qin 
(2003). As shown in Allgower and Zheng (2000), Mayne 
et al. (2000) and Rawlings and Mayne (2009), theory 
for such controllers is also widely understood both for 
linear and nonlinear systems. The control method itself 
deals with the problem of approximately solving an infinite 
horizon optimal control problem which is computationally 
intractable in general. Reasons for its success are on the 
one hand its capability to directly incorporate constraints 
, depending on the states and inputs of the underlying 
' process. On the other hand, the fundamental steps of this 
method are very simple: First, a solution of a finite horizon 
optimal control problem is computed for a given initial 
value. In a second step, the first element of the resulting 
control is implemented at the plant and in the last step, the 
finite horizon is shifted forward in time. As a consequence, 
the method is iteratively applicable and reveals the control 
to be a static state feedback. 

Unfortunately, stability of solution of the infinite hori- 
zon problem may be lost due to considering only finite 
horizons. Over the last two decades, several solutions 
have been proposed to cope with this issue, see, e.g., 
Keerthi and Gilbert (1988), Chen and Allgower (1998) 
and Griine and Rantzer (2008). All these approaches re- 
quire the horizon to be sufficiently long and computing 
the minimal required horizon length is computationally 
demanding. However, the horizon needs to be chosen as a 
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worst case scenario which is usually needed to cope with 
small regions of state space only. Our aim in this work 
is to develop online applicable adaptation strategies for 
the horizon length which guarantee stability of the closed 
loop. In particular, we follow the approach of Griine and 
Pannek (2009) where different suboptimality estimates 
have been developed to measure the performance of the 
model predictive controller. Based on these estimates, we 
propose a simple technique to locally fit the horizon to 
the control task, the current state of the system and also 
to the MPC internal information. Due to the change of 
the structure of the controller, however, known stability 
proofs and suboptimality results cannot be applied. To 
cover these issues, we present a stability result for MPC 
with varying optimization horizons using mild additional 
conditions. To some extend adaptation strategies of the 
horizon are known in the literature, see e.g. Frew et al. 
(2006) and Viquerat et al. (2008), which are heuristics 
based on insight of the specific problem but have shown 
to be applicable in an adaptive model predictive control 
setting. In contrast to that, our approach can be proven 
rigorously and doesnot require any insight into the process 
under consideration (note that different to our intention 
the term adaptive model predictive control is also used to 
incorporate model uncertainties, see, e.g., Michalska and 
Mayne (1993) and Adetola et al. (2009)). 

The paper is organized as follows: In Section 2 we describe 
and motivate the problem setup. Section 3 deals with the 
a posteriori and a priori suboptimalty estimates which will 
be the foundation of our analysis. In the following Section 
4, we show how the stated stability results and estimates 
can be extended to the case of varying optimization hori- 
zons. Thereafter, we state a simple shortening and prolon- 
gation strategy based on the suboptimality estimates given 



in Section 5. In order to show the apphcabihty and effec- 
tiveness of our approach, Section 6 contains a numerical 
example of the adaptive MPC approach. The final Section 
7 concludes the paper and points out directions of future 
research. 

2. SETTING 

Within this work we analyze nonlinear discrete time con- 
trol systems of the form 

x{i + l) = f{x(i),u(i)), x{0)=xo (1) 

with x{i) eX C X and u{i) £ U C C/ for i G No. For the 
considered systems the state space X and the control value 
space U are arbitrary metric spaces. Hence, all following 
results also apply to the discrete time dynamics induced by 
a sampled infinite dimensional system, cf. Ito and Kunisch 
(2002) or Altmiiller et al. (2010). Here, we denote the space 
of control sequences u : No — ^ U by 0^° and the solution 
trajectory for given control u G by a;„(-). Additionally, 
the sets X and U incorporate possible restrictions on the 
state and control respectively. 

In the following, we aim at finding a static state feedback 
u = fi{x) S for a given control system (1) which 
minimizes the infinite horizon cost functional 

oo 

Joo{xo,u) = ^l{Xy,{i),u{i)) (2) 

with stage cost I : X x U ^ Rq . The correspond- 
ing optimal value function is denoted by Voo{xo) = 
inf ^gujNo Joo (xq , u) and throughout this paper we assume 
that the minimum with respect to m S U'^" is attained. 
The optimal value function Voo{-) can be used to define 
the infinite horizon feedback law 

fi{xu{i)) := aTgmm{Voo{xu{i + 1)) + l{xu{i),u)} (3) 

for which one can show optimality using Bellman's op- 
timality principle. Since the computation of the desired 
control law requires the solution of a Hamilton Jacobi- 
Bellman equation, we use a model predictive; control ap- 
proach in order to avoid the problem of solving an infinite 
horizon optimal control problem. The fundamental idea of 
such a model predictive controller is simple and consists 
of three steps which are repeated at every discrete time 
instant during the process run: First, an optimal control 
for the problem on a finite horizon [0, N] is computed given 
the most recent known state of the system xq- Then, the 
first control clement is implemented at the plant and in the 
third step the entire optimal control problem considered 
in the first step is shifted forward in time by one discrete 
time instant which allows for itcrativcly repeating this pro- 
cess. In the literature this method is also termed receding 
horizon, see, e.g., Mayne et al. (2000). 

In contrast to the infinite horizon optimal control (3), the 
problem in the second step is to minimize the truncated 
cost functional on a finite horizon 

Jn{xo,u) = ^ l{xu{k,XQ),u{k)). (4) 

fe=0 

The truncated horizon defines the set of discrete time 
instances I := {0, . . . , TV — 1}. Here, we assume the first 
instant to be denoted by zero for each optimal control 



problem within the MPC problem. In particular, we focus 
on the implementation of a constrained model predic- 
tive controller without additional stabilizing endpoint con- 
straints or a Lyapunov function type endpoint weight, see, 
e.g., Keerthi and Gilbert (1988) and Chen and Allgower 
(1998), respectively. 

Throughout this work, we denote the closed loop solution 
at time instant i by x{i) while Xu{-,xo) denotes the open 
loop trajectory of the prediction. Moreover, we use the 
abbreviations 

un{-, Xq) = argmin Jat (xq, u) (5) 

un{xo) = Un{0,Xo) 

for the minimizing open loop control sequence of the 
reduced cost functional and its first element respectively. 
We call Vn{xo) = min^gpN Jn{xo,u) the optimal value 
function of the finite cost functional (4) and, for notational 
purposes, we use iiAr(i,xo) to represent the i-th control 
value within the open loop control sequence corresponding 
to the initial value xq when it is necessary to distinguish 
between two or more difii'erent open loop controls. Hence, 
if the initial value a;„„ (0, a;o) = a;o is given, then the open 
loop control (5) induces the open loop solution 

Xutiik + l,Xo) = / (x„„(fc,Xo),UAr(fc,.To)) (6) 

for all time instances k on the optimization horizon I \ 
{N}. Similarly to (3), the closed loop control can be 

defined as 

HNix{i)) := argmin{VAf_i(x(i -I- 1)) -I- ^(x(i),w)} (7) 

and the corresponding closed loop system is given by 

x{i + l)=f{xii),f,Nixii))) (8) 

for all i e No. 

Note that due to the truncation of the infinite horizon 

cost functional (2) to the finite MPC cost functional 
(4) , stability and optimality properties of the closed loop 
solution (7), (8) induced by the infinite horizon optimal 
control (3) are not preserved in general. 

Here, our aim is to show that in order to guarantee stability 
of the closed loop (7), (8) for any initial value x € X, 
the requirement of considering the worst case optimization 
horizon N for all initial values a; G X can be weakened. 
Additionally, the resulting closed loop trajectory satisfies 
locally a predefined degree of suboptimality compared to 
the infinite horizon solution (1), (3). 

3. STABILITY FOR STANDARD NMPC 

The measure of suboptimality we consider in the following 
is the difference between the infinite horizon cost induced 
by the MPC law /xjv(-), that is 

DO 

V^-{xo):=J2^i^i')'l^^i^i')))' (9) 

and the finite horizon cost V]\[{-) or the infinite horizon 
optimal value function 14o(')- In. particular, the latter 
give us estimates on the degree of suboptimality of the 
controller /Uiv(-) of the MPC process. For this purpose, we 
make extensive use of the suboptimality estimates derived 
in Griine and Pannek (2009). 



Proposition 3.1. (A posteriori Estimate). Consider a feed- 
back law /ijv : X — !> U and its associated trajectory x{-) 
according to (8) with initial value a:(0) = Xo gX. If there 
exists a function V/v : X Rq satisfying 

VN{x{i)) > VN{x{i + 1)) + al{x{i), (10) 
for some a S (0, 1] and all i € No, then 

aVoo{x{i)) < aV^'^ixii)) < Vn{x{i)) < V^{x{i)) (11) 
holds for all i e Nq. 

Since all values in (10) arc computed throughout the 
NMPC process, a can be easily computed online along the 
closed loop trajectory. Thus, (10) yields a computationally 
feasible and numerically cheap way to estimate the degree 
of suboptimality of the trajectory. 

Due to the fact that VN{x{i + 1)) in (10) is unknown at 
runtime, Proposition (3.1) yields an a posteriori estimator. 
However, we can also utilize a more conservative a priori 
estimate if we assume the following: 

Assumption 3.2. For given N , Nf) e'H, N > Nq> 2, there 
exists a constant 7 > such that for the open loop solution 
{i,x{i)) given by (6) the inequalities 

VNA^u^i.N -No,x{i))) 



7- 



1 



< 



< max l{xu,^{N - j,x{i)),iij-i{xu,^{N - j,x{i)))) 

]=2,...,No 

Vkixu^iN -k,x{i))) 



7 + 1 



< 



< lixu^iN - k,x{i)),Hk(xu^,{N - k,x{i)))) 
hold for all /c e {A^o + 1, • ■ • , N} and all i € Nq. 
Theorem 3.3. (A priori Estimate). Consider 7 > and N, 
No eN, N > No such that (7+1)^-^" > 7^-^0+2 holds. 
If Assumption 3.2 is fulfilled for these 7, A^ and No, then 
the estimate (11) holds for all t G No where 



a := 



(7 + 1) 



N-No 



„iV-Afo + 2 



/ 



(7+1) 



N-Nn 



(12) 



Note that we cannot expect the relaxed Lyapunov in- 
equality (10) or Assumption 3.2 to hold in practice. In 
many cases the discrete time system (1) is obtained from 
a discretization of a continuous time system, e.g. sam- 
pling with zero order hold, see Nesic and Teel (2004). 
Hence, even if the continuous time system is stabilizable 
to a setpoint x* and no numerical errors occur during 
optimization and integration, the corresponding sampled- 
data system is most likely practically stabilizable at x* 
only. However, suboptimality results can be extended to 
cover the case of practical stability as well, see Griine 
and Rantzer (2008) and Griine and Pannek (2009). Since 
extending the stability results we will present now to cover 
the practical case can be done analogously, see Pannek 
(2009), we restrict ourselves to the case of asymptotic 
stability for simplicity of exposition. 

4. STABILITY UNDER ADAPTATION 

As stated at the end of Section 2, wc aim at weakening the 
worst case nature of the optimization horizon N. Here, one 
has to keep in mind that if a model predictive controller 
shall be designed for a given application, then stability 
of the resulting closed loop (8) needs to be guaranteed 



for the entire working range X. In practice, this may lead 
to very large optimization horizons N. Yet, most points 
visited by the closed loop (8) wc do not require such a 
large optimization horizon in order to guarantee stability. 

Here, wc focus on locally guaranteeing a decrease of the 
cost function for each step of the MPC process and modify 
the horizon length N to fulfill this task. Similar to the 
suboptimality results from Section 3, we want to measure 
this decrease in terms of the running cost ?(•, •) such that 
a given suboptimality bound a G (0, 1) is locally satisfied. 

Since we are now dealing with varying optimization hori- 
zons, we intuitively extend our notation from Section 3 
by adding the used optimization horizon as an argimient, 
i.e. a{N) denotes the suboptimality degree a with horizon 
N. Moreover, since the resulting closed loop control now 
depends on a sequence {Ni)i^^ we denote such a control 
law by /i(jVi)- 

An abstract adaptive MPC algorithm which locally ac- 
complishes the task of guaranteeing a decrease in the cost 
function is the following: 



(1) Given x{i) and Ni do 

(la) Compute optimal control on horizon Ni 

(lb) Compute suboptimality degree a{Ni) 

(Ic) If a(iVj) > a: Call shortening strategy for Ni 

Else: Call prolongation strategy for Ni 
while a{Ni) < a 

(2) Implement the first control component iipf.{x{i)) := 
u{0, x(i)) 

(3) Set i := i + 1 and shift the optimization horizon 
forward in time 



In this context, we distinguish the following degrees of 
suboptimality: 

Definition J^.l. (Suboptimality Degree), (i) Given a set X, 
then we call a := max{a | (10) holds Vx(n) = a; € X} the 
global suhoptimality degree. 

ill) Given a point a; G X, then we call a := max{a | 
(10) holds for x(n) = x] the local suhoptimality degree. 
(iii) Given a closed loop trajectory x{-) we call a := 
max{a | (10) holds Vn e No} the closed loop suhoptimality 
degree. 

The problem which we are facing for such an adaptive 

MPC algorithm is the fact that none of the existing sta- 
bility proofs, see, e.g., Keerthi and Gilbert (1988), Chen 
and Allgower (1998), Grimm et al. (2005), Jadbabaie and 

Hauser (2005), Griine and Pannek (2009) and Griine et al. 
(2010), can be applied in this context since these results 
assume N to be constant while here the optimization hori- 
zon Ni may change in every step of the MPC algorithm. 

The major obstacle to apply the idea of Proposition 3.1 in 
the context of varying optimization horizons N is the lack 
of a common Lyapunov function along the closed loop. To 
compensate for this deficiency, we make the following mild 
assumption: 

Assumption 4-2. Given an initial value x € X and a 
horizon length N < 00 such that /xjv(-) guarantees local 
suboptimality degree a{N) > a, a G (0,1), we assume 
that ioT N > N, N < 00, there exist constants Ci,Ca > 



such that the inequaUties 
CiI{x,ij,n{x)) < l{x,fj,~{x)) 



V~{x)-V~inx,^^Nix))) 



where we defined 



N 



V^ix)-V~{f{x,f,~{x)) 



N 



hold where a{N) is the local suboptimality degree of the 
controller Atj^(-) corresponding to the horizon length N. 

Note that Assumption 4.2 is indeed very weak since 
for one we allow for non -monotone developments of the 
suboptimality degree a(-) if the horizon length is increased 
which may occur as shown in Di Palma and Magni (2007). 
Here, we only make sure that if a certain suboptimality 
degree a € (0, 1) holds for a horizon length A'', then the 
estimate a{N) does not drop below zero if the horizon 
length N is increased. 

Considering the value of l{x, ii~{x)), we notice that it may 

tend to zero if N is increased, hence we have that C; is 
in general unbounded. The special case 1{x,ijl~{x)) = 0, 
however, states that the equilibrium of our problem has 
been reached and can be neglected in this context since 
this implies 1{x,Hn{x)) =0 allowing for arbitrary Ci. 

Given Assumption 4.2, we obtain stability and a per- 
formance estimate of the closed loop in the context of 
changing horizon lengths similar to Proposition 3.1. 

Theorem 4.3. (Stability of Adaptive MFC). Consider a e 
(0,1) and a sequence {Ni)i^^g, Ni € N, where N* = 
max{A^j I i E N}, such that the MFC feedback law /^(at.) 
defining the closed loop solution (8) guarantees 

VnM^)) > VnAx{^ + l)) + alixii),tXNMi))) (13) 

for all i € No- If additionally Assumption 4.2 is satisfied 
for all pairs of initial values and horizons {x{i), N^), z G Nq, 
then we obtain 

acVoo{x{n)) < acV^^^'^xin)) < Vjv(a;(n)) < Voo{x{n)) 

to hold for all n € No where ac '■= min a. 

jeN>„ ' 

Proof: Given a pair {x{i), Ni), Assumption 4.2 guarantees 

a{Ni) < a{N)/cL'^ for N > Ni. Now we choose N = 
N* within this local suboptimality estimation. Hence, 

we obtain a < a{Ni) < a{N*)/Ca'' using the relaxed 
Lyapunov inequality (13). Multiplying by the stage cost 
l{x{i), fiNi{x{i))), we can conclude 

al{x{i), ^N.{x{i))) < 
a{N*). 



< 



a 



-l{x{i),iiNi{x{i))) 



< 



VN*{x{i))-VN*{f{x,tMN*{x))) 

&::h{x{i),tiN*{x(i))) 

VN*{x{i)) - VN*{f{x,HN,{x))) 
Mi) Mi) 



l{x{i),iiNi{x{i))) 



using (13) and Asstunption 4.2. Summing the running 
costs along the closed loop trajectory reveals 



ac 



K 



l{x{i), ix{i))) < Vn* ix{n)) - Vn* {x{K + 1)) 



Since Vn* {x{K+1)) > holds, taking K to infinity reveals 



K 



acK^*"''(a;("-)) = ac lim 'S^l{x{i), iJLNi{x{i))) 

K^oo ^ — ' 



< VN*{x{n)) 

Since the aVoc{x{n)) < acVx'^'\x{n)) and VN*{x{n)) < 
Voo{x{n)) hold by the principle of optimality, the assertion 
follows. □ 

Comparing Proposition 3.1 and Theorem 4.3, we see that 
the closed loop estimate ac may be smaller than the 
local suboptimality bound a but due to C; , Cq, > we 
can guarantee ac > 0. Yet, ac may become very small 
depending on Ca and C; from Assumption 4.2. 

5. A SIMFLE ADAFTATION STRATEGY 

Since now we have shown asymptotic stability of a MFC 
closed loop trajectory with varying optimization horizon, 
we show a very simple approach to guarantee the local 
suboptimality requirement a{Ni) > a. To this end, we 
assume the system to be controllable, i.e. 
Assumption 5.1. Given a e (0,1), for all € X there 
exists a finite horizon length N = N{xo) & ^ such that the 
relaxed Lyapunov inequality (10) holds with a{N) > a. 

Theorem 5.2. (Shortening Strategy). Consider an optimal 
control problem (5), (6) with initial value x(«), horizon 
A'i £ N and fixed suboptimality bound 75 G (0,1) and 
denote the optimal control sequence by u* . Suppose there 
exists an integer i G No, < i < A^i such that 

{k,x{i))) - VN,-k(xu*{k + l,x{i))) > 
> al{xu* (k, x{i)), UNi-kixu* [k, x{i)))) (14) 

holds true for all < A: < i. Then, setting iV^+fe = Ni — k 
and fj,Ni+k{^{i + k)) = u*{k) for < A; < i — 1, inequality 
(13) holds for fc = i, . . . , i + i — 1 with a = a. 

Proof: The proof follows directly from the fact that for 
A'Afi+fc(a;(« + fc)) = u*{k) the closed loop trajectory satisfies 
x{i + k) = Xu*{k,x{i)). Hence, (13) follows from (14). 

With the choice A^i+fe = Ni — k, due to the principle of 
optimality we obtain that the optimal control problems 
within the next i — 1 NMFC iterations are already solved 
since niq.^^{x{i + k)) can be obtained from the optimal 
control sequence w*(-) G i]'^ {x{i)) computed at time i. 
This implies that the most efficient way for the reducing 
strategy is not to reduce Ni itself but rather to reduce the 
horizons A^i+fc by k for the subsequent sampling instants 
i-|-l,...,i+Z;i.c.,wc choose the initial guess of the horizon 
A^j+i = A^j — 1. Still, if the a posteriori estimate is used, 
the evaluation of (14) requires the solution of an additional 
optimal control problem in each step. 
In order to to use the a priori estimate given by Theorem 
3.3 the following result can be used as a shortening 
strategy: 

Theorem 5.3. (A priori Shortening Strategy). Consider a 
optimal control problem (5), (6) with initial value x{i) 
and horizon Ni, N G N, Ni > N > 2 and denote the 
optimal control sequence by u*. Moreover, the suboptimal- 
ity bound a G (0, 1) is fixed inducing some 7(-) via (12). 



Suppose there exists an integer i Q No, < i < Ni — Nq — I 
such that for all < k < i there exist ji < j{Ni — k) 
satisfying 



VNAxu*{Ni-N,x(i))) 



< 



l^ + 1 

< max l{xu*{Ni- j,x{i)),iJ,j-i{xu*{Ni- j,x{i)))) 

j=2,...,N 

Vk,{xu* {N,-ki,x{i))) 



(15) 



< 



7i + 1 

< l{xu* {Ni - ki,x{i)), {xu* {N -ki, x{i)))) 



(16) 



for all ki e {iV+1, . . . , Ni-k}. Then, setting Ni+k = Ni-k 
and ^j\j.^^{x{i + k)) = u*{k) for < fc < i — 1, inequality 
(13) holds for A; = 2, ...,« + i — 1 with a = a. 

Proof: Since (15), (16) hold for k = 0, Theorem 3.3 
guarantees that the local suboptimality degree is at least 
as large as a. If i > holds, we can make use of the fact 
that for /ZAr.^j^(x(i + fc)) = u*{k) the closed loop trajectory 
satisfies x{i + k) = x„*(fc,x(i)). By (15), (16), we obtain 
Assumption 4.2 to hold for k = i, . . . Accordingly, 
the assertion follows from Theorem 3.3 which concludes 
the proof. □ 

Note that while the a priori estimate from Theorem 3.3 is 
slightly more conservative than the result from Proposition 
3.1, it is also computationally less demanding if the value 
Nq is small. 

In contrast to this efficient and simple shortening strategy, 
it is quite difficult to obtain efficient methods for prolon- 
gating the optimization horizon iVj. In order to obtain a 
simple prolongating strategy, we invert the approach of 
Theorem 5.2, i.e. we iteratively increase the parameter A'^ 
until the requirement a{Ni) > a is satisfied. 

Theorem 5.4- (Prolongation Strategy). Consider an opti- 
mal control problem (5), (6) with initial value x{i) and 
Ni e N. Moreover, for fixed a G (0, 1) suppose Assump- 
tion 5.1 to hold. Then, any algorithm which iteratively 
increases the optimization horizon iVj terminates in finite 
time and computes a horizon length Ni such that (13) 
holds with local suboptimality degree a. 

Proof: Follows directly from Assumption 5.1. □ 

Unfortunately, if (13) docs not hold, it is in general difficult 
to assess by how much N should be increased such that 
(13) holds for the increased Ni. The most simple strategy 
of increasing Ni by one in each iteration shows satisfactory 
results in practice, however, when starting the iteration 
with Ni, in the worst case it requires us to check (13) 
N — Nn + 1 times at each sampling instant. In contrast to 
the shortening strategy, the principle of optimality cannot 
be used here to establish a relation between the optimal 
control problems for different Ni and. moreover, these 
problems may exhibit different sohition striictures which 
makes it a hard task to provide a suitable initial guess for 
the optimization algorithm. 



6. NUMERICAL RESULTS 

To illustrate the eff'ects of using an adaptive NMPC, we 
consider a highrack warehouse 



v{t)=U2it) (17) 

Mt) = -kkt) - -(t)sm{6it)) -uiit) cosUU)) 

V 

where for simplicity of exposition the rope is modeled 
as a pendulum with variable length. Here, x denotes the 
position of the crab along the highrack, v represents the 
length of the rope of the crane and (p corresponds to the 
angle of deflection of the rope. Moreover, g = 9.81 and 
fc = 0.1 denote the gravitational constant and the inertia 
of the angle of the rope, respectively. 

For this example, we use MPC to generate a feedback for a 
representative transport action of a pallet from xo = —3, 
uo = 5 to Xref = 3, Vref = 2 (with zcro derivatives in 
initial and target position) while maintaining the state and 
control constraints X = [-5, 5]^ x [1, 4] x [-1, 2] x [-1, l]xR 
and U = [—5, 5] x [—1, 2]. To this end, we use the running 
cost 

T 

l{x,u) = j ci4>'^{t)v'^{t) + C2gv{t){l - cos(0(i))) 

+ caixit) - X)' + C4X^{t) + C5{v{t) - vf 

+ CQi?{t) + c-j {vl{f) + ul(t)) dt 

with constants c\ = 0.25, C2 = 0.5, C3 = 40, 04 = 05 = 
C6 = 20 and C7 = 0.1 and the sampling period T = 0.2. To 
solve the optimal control problem arising throughout the 
MPC procedure, we use a direct approach, i.e. discretize 
the continuous time problem and use an SQP method to 
solve the resulting optimization problem. Here, we set the 
tolerance levels toloDE = 10~^ and tolggp = 10~® for the 
differential equation solver and the optimization method 
respectively. 

Since the adaptive MPC algorithm allows us to set the 
lower bound of the degree of suboptimality a directly, we 
first investigate the a-depending quality of a controlse- 
quence on the closed loop cost V^"(a;o). To this end, we 
terminate the algorithm when the condition l{x{t), u{t)) < 
10~^ is satisfied. The data we obtained for this setting is 
displayed in Figure 1. 

y^^(a;o)-lo3 
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Fig. 1. Development of V^" {xo) for different suboptimality 
bounds a 

Here, one can nicely observe that the closed loop costs 
caused by the adaptive MPC feedback /U(jVi) are decreasing 
as the lower bound a is enlarged. This is right the behav- 
ior one would expect from the theoretical construction. 



However, using the adaptive MPC approach, larger a- 
values not only provides a much better control sequence in 
terms of generated costs. Wo also like to mention that the 
total simulation time required to satisfy the termination 
criterion is also decreasing as a is enlarged which is due to 
the use of larger optimization horizons Ni throughout the 
run of the simulation. 

In Figure 2, we additionally plotted the optimization 
horizon sequences {Ni) for the selected values of a. This 
figure demonstrates clearly the horizon incrementations 
during acceleration- and deceleration-phases. In partic- 
ular, a large optimization horizon is required to satisfy the 
desired decrease in the relaxed Lyapunov inequality (13) 
upon start of the simulation run which is then reduced as 
the crab moves towards its destination. In order to reduce 
the possibly occuring overshoot, the method automatically 
increases the horizon again. During the final leveling phase, 
again no large horizons are needed to satisfy (13). 
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Fig. 2. Horizon length N during simulation with a = 
0.6(solid) and a = 0.2(dashed) 

In Figure 2, one can also see that the spike in the horizon 
length occurs earlier for a = 0.6. Again, this corresponds 
to the MPC procedure recognizing the possible overshoot 
by means of (13). 

7. CONCLUSION 

In this work we have shown stability and suboptimality 
estimates for model predictive controllers with varying 
optimization horizon. This result allows for developing 
strategies to adapt the horizon length instead of using a 
worst case estimate and control the quality of the feedback 
directly. 

Future work may concern reducing the computational 

effort required to evaluate the suboptimality estimates. 
Moreover, development and investigation of alternatives 
to prolongate the optimization horizon will be an issue, 
i.e. by combining information of several iterates. 
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